Hydrodynamics of the Kuramoto-Sivashinsky Equation in Two Dimensions 
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The large scale properties of spatiotemporal chaos in the 2d Kuramoto-Sivashinsky equation 
are studied using an explicit coarse graining scheme. A set of intermediate equations are obtained. 
They describe interactions between the small scale (e.g., cellular) structures and the hydrodynamic 
degrees of freedom. Possible forms of the effective large scale hydrodynamics are constructed and 
examined. Although a number of different universality classes are allowed by symmetry, numerical 
results support the simplest scenario, that being the KPZ universality class. 
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A major goal in the study of spatiotemporal chaos 
(STC) jl| is to obtain quantitative connections between 
the chaotic dynamics of a system at small scales and 
the apparent stochastic behavior at large scales. The 
Kuramoto-Sivashinsky (KS) equation S 



9th 



-V 2 h-V A h+{Vh) 2 
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has been used as a paradigm in efforts to elucidate the 
micro- macro connections |3|,pJ . 

The qualitative behavior of the KS equation is quite 
simple. Cellular structures are generated at scales of the 
order £q — 2\f2-K due to the linear instability. These cells 
then interact chaotically with each other via the non- 
linear spatial coupling to form the STC steady state at 
scales much larger than £q. The characterization of the 
STC state has been studied extensively in one spatial 
dimension It was conjectured early by Yakhot ||, 

based partially on symmetry grounds, that the large scale 
behavior of the one dimensional KS (ld-KS) equation 
is equivalent to that of the Id noisy Burgers equation, 
also known as the KPZ equation H. This conjecture has 
since been validated by detailed numerical studies 
More recently, an explicit coarse-graining procedure was 
used by Chow and Hwa M to derive a set of coupled ef- 
fective equations describing the interaction between the 
chaotic cellular dynamics and the long wavelength fluc- 
tuations of the /i-field. From this description, the large 
scale (KPZ-like) behavior of the ld-KS system can be 
predicted quantitatively from the knowledge of various 
response functions at the 'mesoscopic scale' of several 
Vs. 

The behavior of the 2d-KS equation is not as well un- 
derstood. The simplest scenario is the generalization of 
Yakhot's conjecture to 2d, with the large scale behavior 
described by the 2d-KPZ equation, 



d t h = vV 2 h+^(Vh) 2 +r l (r,t), 



(2) 



the asymptotic scaling properties of Eq. (|2j) are described 
by 'strong-coupling' behavior with algebraic (rather than 
logarithmic) scaling in the roughness of h and super- 
diffusive dynamics The length scale at which the 
asymptotic regime is reached is given by £ x ~ e 87T ^ 9 , 
where g = X 2 D/u 3 . At scales below £ x , the effect of the 
nonlinear term in (||) can be accounted for adequately via 
perturbation theory. The system behaves in this 'weak- 
coupling' regime as a linear stochastic diffusion equation 
with additive logarithmic corrections 

Previous studies of the 2d-KS equation |l(]-[l^] found 
behavior consistent with linear diffusion with loga- 
rithmic corrections but had different interpretations. 



where v > can be interpreted as a stabilizing 'surface 
tension' for the height profile h, and 7/ a stochastic noise 
with (r/{r,t)ri(r',t')) = 2DS 2 {r - r')S(t - *')■ For v > 0, 



Jayaprakash et al. |12| performed a numerical analy- 
sis akin to Zaleski's on the ld-KS Q, and concluded 
that their results were consistent with the weak-coupling 
regime of the 2d-KPZ equation, with (in principle) a 
crossover to strong-coupling beyond a length of £ x s» 
10 26 £o, for 9 = 0.4. Procaccia et al. 0Jll| used a com- 
parative Dyson- Wyld diagrammatic analysis of the two 
equations to argue that 2d-KS and 2d-KPZ cannot be- 
long to the same universality class |]l3| . They maintained 
instead that the asymptotic behavior of the 2d-KS equa- 
tion is described by a 'nonlocal' solution, consisting of 
diffusion with multiplicative logarithmic corrections pi] . 
We feel that the ensuing debate jll| failed to rule out 
either interpretation. 

It is very difficult to distinguish between the above 
two scenarios numerically, as one must resolve different 
forms of logarithmic corrections to the (already logarith- 
mic) correlation function of the linear diffusion equation. 
Theoretically, there is no a priori reason why simple sym- 
metry considerations such as Yakhot's should be valid in 
two and higher dimensions. Unlike in Id where there are 
only scalar density fluctuations, the 2d case is compli- 
cated because three or more large-fc modes can couple 
and contribute to low- A: fluctuations. Such nonlocal in- 
teractions in k may not be adequately accounted for in 
the type of analysis performed in Refs. ^,^2|, which nu- 
merically impose KPZ dynamics and then test for self- 
consistency. 
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In this study, we perform a systematic symmetry anal- 
ysis, taking into account the possibility of large-k cou- 
pling. Specifically, we extend the coarse graining proce- 
dure of Ref. |6| to two dimensions to derive a set of cou- 
pled equations describing the local arrangement of cells, 
and study their effect on the macroscopic dynamics of 
the ft,-field. The resulting behavior depends crucially on 
the small scale arrangement of the cells. In the sim- 
plest case, the strong-coupling 2d-KPZ behavior is re- 
covered. Nevertheless, more complicated behaviors are 
allowed if the microscopic cellular arrangement exhibits 
spontaneous rotational symmetry breaking. A number of 
possible scenarios are listed for this case. To determine 
which of the allowed scenarios is selected by the 2d-KS 
equation, we performed numerical measurements of the 
cellular dynamics at the mesoscopic scale of 4 to 16 £q's. 
Our results disfavor the occurrence of the more exotic 
scenarios, leaving the strong-coupling 2d-KPZ behavior 
as the most likely possibility. 

As in Id, we coarse grain over a region of size L x L, 
where L is several times the typical cellular size Iq. h(r, t) 
is separated into fast cellular modes h > and slow long 
wavelength modes h < . Inserting h(r,t) = h < (r,t) + 
h > (r,t) into Eq. (|l|), we obtain the following equations 
for the fast and slow modes: 



d t h>= -V 2 /i> 



V 4 /i> -I- (V/i>)| + 2{Wh > ■ V/i<) (3) 



dth<= -S7 2 h < + (V/i<) 2 + w(r, t) + 0(V 4 /i<). 



(4) 



where w(r, t) = (V/i>) 2 , is the only contribution of the 
fast modes on the dynamics of h < . It can be interpreted 
as the 'drift rate' of h < over a regime of L x L centered at 
r. To specify the dynamics of h < , it is necessary to ob- 
tain the dynamics of w from the fast mode equation (|^) . 
Due to the structure of the nonlinear term, we must con- 
sider the tensor W, with elements Wij — 2 dihy ■ djhy 
where the over-line denotes a spatial average over the 
coarse-graining scale L. It is convenient to introduce the 
curvature tensor K, with elements Kij — 2didjh < . In this 
notation, w = ±TrW and n = V 2 /i< = ^Tr K. Taking 



the time derivative of W%j and using (|3j) , we obtain 
<9 t W = F[W]+W-K + K- W 



(5) 



where F [W] contains purely fast mode dynamics and will 
be described shortly. The forms of the last two terms in 
Eq. (H) are fixed by the Galilean invariance of the KS 
equation and are exact. 

Equation (||) can be made more transparent by rewrit- 
ing the two tensors as W = w ■ 1 + w ■ Q(4>) and 
K = k ■ 1 + K ■ Q(0), where 1 is the identity ma- 
trix and Q(a) is a unit traceless matrix, represented 
by an angle a, e.g., Q\2{a) = Q2i(ot) = sin(2a) and 
Qn(a) = —Q22(a) = cos(2a). Adopting vector nota- 
tion if) = (w cos 2(f), w sin 20) and x = (k cos 29, k sin 29), 
Eq. (El) can be rewritten as 



dtw = f[w] + 2kw + 2\ ■ if) 
d t if> = ip + Kif) + wx, 



(G) 
(7) 



to leading order, with / and tp obtained from the appro- 
priate decomposition of F. 

Eqs. @ and (Q), together with the slow mode equation 
(Q), form a closed set of coarse grained equations which 
specifies the dynamics of h < once the effective forms of 
the small scale dynamics, i.e., / and ip are given. These 
equations are constructed from symmetry considerations, 
and can be regarded as the more complete generalization 
of Yakhot's conjecture for two dimensions. We first dis- 
cuss the physical meaning of the coarse grained variables 
appearing in W and K. 

The tensor K describes the local curvature of the slow 
modes h < . With k = 0, we have a symmetric paraboloid 
- a 'valley' if n > or a 'hill' if k < 0. With k = 0, 
we have a 'saddle', with k and 9 specifying the strength 
and the orientation. The tensor W characterizes the local 
packing of the cells. As in the Id case [0 , w gives the lo- 
cal cell density. The traceless component of W describes 
the local anisotropy in cell packing. Fig. [l](a), shows an 
example of an arrangement of cells where ip is pointed 
in the y-direction. Fluctuations in the anisotropic part 
of the curvature K will affect the local cell packing. For 
example, the cell density at the bottom of a valley will 
be higher and a saddle configuration in h < will induce 
anisotropy. Eqs. (^|) and (0) describe these effects of 
curvature quantitatively, much like the relation between 
stress and strain in elastic systems. Cell packing in turn 
influences the slow mode dynamics via the w term in 
Eq. (^). The anisotropic parts of K and W are invariant 
upon a rotation by 180°; see Fig. 0(a). Thus, we can 
view the vector field if)(r, t) as a 'nematic' order parame- 
ter describing the local cellular orientation, and x as an 
applied field biasing if towards a specific orientation. 

If we turn off the applied field k and x i n Eqs. (p) and 



for each coarse- 



(0), we have d t w = f[w] and d t ip = <p 
grained region; thus / and <p describe the small scale dy- 
namics. Even for a coarse-grained region of a few ^o's, 
the small scale dynamics of h are already chaotic. The 
fields w(t) and if>(t) are 'projections' of this small scale 
chaotic dynamics. They can be quantitatively character- 
ized numerically as we will present shortly. Before doing 
so, we first construct some possible scenarios. 

We expect that the ft,-field has on average a finite 
drift rate, i.e., a finite time-averaged value of w. The 
simplest dynamics of w is then dtw(r,t) — f[w] = 
— a(w — wo) +£(r, t), where £ is a stochastic forcing mim- 
icking the chaotic small scale dynamics, and wq is a con- 
stant. This yields w(r,t — * 00) — * wq. The behavior of 
if) is less straightforward. In the simplest scenario, we 
have dtif) = —aip + (f(r, t) to leading order, with £ being 
a vectorial stochastic forcing. Eq. (]?]) then yields (in the 
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hydrodynamic limit) ip ~ (2wo/a)x, where we took the 
asymptotic result w = wq and assumed that the typi- 
cal curvature k is small. Note that in this scenario, the 
cellular orientation passively follows the curvature. In 
particular, there is no oricntational anisotropy on aver- 
age if there is no external forcing. Inserting this result 
and f[w] into (^|), we find in the hydrodynamic limit 



w ~ Wo + ^V 2 h < + ^ + 0((d l d 1 h < ) 2 ). (8) 
a 



Substituting (||) into (0) yields an equation for h < of the 
KPZ form (j^) to leading order, with v = (2wq/cx) — 1 
and r)(r,t) = £(r,t)/a. Dynamics of the KPZ univer- 
sality class will be obtained if v > and the noise £ is 
uncorrelated between different coarse-graining regions. 

Unlike the constant a however, there is no a priori rea- 
son why the constant a cannot be negative. This would 
be the case if the microscopic chaotic dynamics has a 
preference for the spontaneous breaking of local isotropy. 
If a < 0, then the dynamics of ip would be more compli- 
cated. Higher order terms, e.g., |^>| 2 ^>, will be needed for 
stability. The minimal equation for (Q) becomes 

d t tp = -aip - (3 ip ip + lV 2 ip + w X + ((r,t) (9) 

where @ is a positive constant, and the 7 term de- 
scribes the coupling of neighboring coarse grained re- 
gions. Eq. (^) describes the relaxational dynamics of a 
nematic liquid crystal under an applied 'field' x- Its be- 
havior depends crucially on the dynamics of the phase 
field, 4>, which is the Goldstone mode associated with 
symmetry breaking. The latter in turn depends on the 
parameters of Eq. (^|), particularly the coupling constant 
7 and the amplitude of the noise £• The possibilities 
along with the effects on h < are: 

Case i) If the noise £ dominates over the spatial cou- 
pling 7, then the local anisotropy will be destroyed at 
large scales due to the proliferation of topological defects 
(disclinations) in <j>. Isotropy is restored and the KPZ 
universality class is recovered. 

Case ii) If the spatial coupling is large, then the di- 
rection of ip may 'phase-lock' with the direction of x, as 
manifested by ((9 — <p) 2 ) <C 1. Solving for the steady 

state of w in this case gives w ~ wq + —k + leading 
to a slow mode equation which is explicitly not KPZ-like 

1 /2 

since k = [(h xx — h yy ) 2 + 4h xy \ ■ (In the KPZ case, k 
comes in at second order and is presumed irrelevant; see 
Eq. ®). 

Case iii) For intermediate parameters, there may ex- 
ist a 'spin wave' phase characterized by (<p(r)<p(0)) — 
alog|r|. Here, 'spin wave' fluctuations would add a 
long range component to the effective KPZ noise, since 
(cos(0(r) — 0(0))) ~ r~ a . For sufficiently small a, it 
would yield dynamics that are not in the KPZ universal- 
ity class. 



To distinguish between these scenarios, we numerical 
measured the 'fast mode' dynamics with simulations on 
systems of size Lx L with L ranging from 32 to 128 with 
periodic boundary conditions. (A single cell had a length 
la « 8.9.) We used a simple spatial discretization scheme 
with an Euler time step of .02. We first checked for 
spontaneous symmetry breaking by examining the dis- 
tribution of ip. We measured ip spatially averaged over a 
system of size L — 32, sampling at time intervals of t = 4 
over a total period of t = 4.8 x 10 5 . The distribution 
of \ip\ = w (normalized by the cylindrical area) is shown 
in Fig. |I](b). The Gaussian shape (within ~ 3%), indi- 
cates a lack of spontaneous symmetry breaking strongly 
supporting the KPZ scenario. 
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FIG. 1. a) Example of a cellular arrangement where ip is 
nonzero with orientation 2(f> = 90°, reflecting a larger 'com- 
pression' in the y direction, b) Probability density of w = \ip\. 
Solid line is a Gaussian fit. 

We next measured the response of W to an imposed 
curvature tensor K using methods similar to those de- 
scribed in Ref. ||. In an STC steady state, we abruptly 
turned on a forcing term in Eq. (|^) of the form h < {x, y) — 
ck~ 2 smk x xsmk y y, with k x = k y = tt/2L, for a range 
of amplitudes c. This configuration yielded eight sepa- 
rate overlapping regions of size L/2 x L/2: four hill and 
valley regions, and four saddle regions with orientations 
9 = 7r/4 and 9 = 37r/4. To account for the spatial de- 
pendence of the curvature K within each forcing region, 
we simply took the imposed curvature to be the average 
curvature. This gave k = V 2 /i< = (4/7r 2 )c. We zeroed 
the first three Fourier modes of /i> after each time step to 
remove the nonlinear slow mode response to the forcing 
(see Ref. [§). 

We found that w only responded to k and ip only 
responded to x, further validating the KPZ scenario. 
Fig. ||(a) shows an example of the time dependent re- 
sponse of the drift rate w to the forcing. The saturated 
amplitude Al(c) was a linear function of c as seen in 
Fig ^|(b). In accordance with Eq. (Q), the ratio 2wo/a 
was computed from the slope, from which we obtained 
v = 14.9 ±0.5 for L — 64 and v = 14.5 ±0.5 for L = 128. 
The averaged response times were oT x = 3.3 ± 0.2 for 
L = 64, and a" 1 = 4 ± 0.2 for L = 128. Similar behav- 
iors were obtained in the hyperbolic forcing region. We 
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verified that ifi oc x, an d found that a 1 = 2.0 ± 0.1 for 
L = 64, and 5" 1 = 3.0 ± 0.3 for L = 128. 
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it belongs to the KPZ universality class. However, for 
all practical purposes, the behavior is well described by 
a stochastic diffusion equation with logarithmic correc- 
tions. 

This research is supported in part by NIH grant KOI 
MH1058 (CC), ONR grant N00014-95-1-1002 (TH), and 
the A.P. Sloan Foundation (CC,TH). 



FIG. 2. (a) Response of the drift velocity w to forcing with 
c = 0.1 applied at t = for L — 128. (b) Saturated amplitude 
Al(c); the slope is 3.15 ± 0.06. 

To characterize the effective stochasticity, we measured 
the two-point correlation function Cl(t) = ((w(t + r) — 
w(t)) 2 ) of the drift rate w for systems of sizes ranging 
from L = 32 to L — 128. The relevant quantity is the 
correlator D l {t) = [C L (oo) - C L (r)]/2 (see Fig. |(a)). 
For short-range correlated noise, the effective noise am- 
plitude D — L 2 f^° Dl{t)cIt is expected to be indepen- 
dent of L. Fig. |3(b) shows measurements of D, with the 
average D — 89 ±5. The numerical values of the effective 
parameters D and v extracted using our coarse-graining 
scheme are in reasonable agreement with that of Ref. [fL2| . 
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FIG. 3. (a) Correlator for L 
different L's. 
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32; (b) Noise amplitude for 



In summary, we have performed an analysis of the 2d- 
KS equation at the mesoscopic scale of several cell sizes. 
By using an explicit coarse graining scheme, we con- 
structed various possibilities for the effective equation of 
motion for the slow modes h < . We found, as pointed out 
before in |l(],[ll[], that the KPZ universality class is not 
the only possibility for the 2d-KS equation. Since our 
analysis is confined to the mesoscopic (large-k) limit, the 
various scenarios obtained are all nonpcrturbative and 
nonlocal in k. The more interesting scenarios involve 
spontaneous breaking of rotational symmetry. Although 
we found no symmetry breaking for the KS equation nu- 
merically, our analysis indicates that such solutions are 
allowed by symmetry, and may occur in other KS-like 
system. For the KS equation proper, we conclude that 
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